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Recent experiment [S.I. Dorozhkin et al., Phys. Rev. Lett. 102, 036602 (2009)] on quantum 
Hall structures with strongly asymmetric contact configuration discovered microwave-induced pho- 
tocurrent and photovoltage magnetooscillations in the absence of dc driving. We show that in an 
irradiated sample the Landau quantization leads to violation of the Einstein relation between the dc 
conductivity and diffusion coefficient. Then, in the presence of a built-in electric field in a sample, 
the microwave illumination causes photo-galvanic signals which oscillate as a function of magnetic 
field with the period determined by the ratio of the microwave frequency to the cyclotron frequency, 
as observed in the experiment. 
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I. INTRODUCTION 

Recent developments in the theory of nonequilibrium 
magnetotransport of a two-dimensional electron gas 
(2DEG) in high Landau levels are motivated by the dis- 
covery of several novel kinds of quantum magnetooscilla- 
tions induced by microwave radiation,—"^ by strong di- 
rect current r^i^i^^i^i^i or by phonons.22i^ii£iiL A par- 
ticular attention has been attracted by the microwave- 
induced resistance oscillations (MIRO)ii^ governed by 
the ratio lo /lUc of the circular radiation frequency w and 
the cyclotron frequency cjc — \e\B/mc (here B is the 
magnetic field and m the effective electron mass). Fur- 
ther experiments on MIRO led to spectacular observation 
of the "zero resistance states" (ZRS) in which the dissi- 
pative components of both the resistance and conductiv- 
ity tend to zeroi^i^i^i^ii These states were explained in 
Ref. [13 as a result of instability leading to formation of 
domains carrying non-dissipative Hall current. 

Initially MIRO were attributed to the displacement 
mechanism which accounts for spatial displacements of 
semiclassical electron orbits due to radiation-assisted 
scattering off disorder,— i^i^i^ Due to Landau quanti- 
zation leading to periodic modulation in the density of 
states (DOS) v{e) ~ v{£ + u)c), the preferred direction of 
such displacements with respect to symmetry-breaking 
dc field oscillates with w/wc- This results in MIRO with 
the phase and period observed in Refs. [l|@HllH@. Later 
it was realized that the dominant contribution to MIRO 
in Refs. T,'2"3"'?,5,'6' is due to ^Hnelastic" mechanism asso- 
ciated with radiation-induced changes in occupation of 
electron states f^i^iiiSiiiS while the displacement mecha- 
nism can be relevant at higher temperatures and only if 
sufficient amount of short-range impurities is present in 
the system^i^i^ or else, at a very strong dc field''^ or 
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microwave power i-s. 

So far, the theoretical research on nonequilibrium mag- 
netooscillations in high Landau levels has been concen- 
trated on the properties of systems which are spatially 



homogeneous on the macroscopic scale42r— Here we de- 
velop more general transport theory applicable also for 
nonuniform carrier and field distributions. From the ex- 
perimental side, the present study is motivated by recent 
experiment^'* which discovered alternating-sign magne- 
tooscillations of photocurrent and photovoltage induced 
by microwaves in the absence of dc driving. The magne- 
tooscillations with a phase and period similar to MIRO 
were observed in a 2DEG with a strongly asymmetric 
contact configuration. The effect was related to the ex- 
istence of built-in electric fields in a sample in thermo- 
dynamic equilibrium, in particular, in vicinity of doped 
contacts. As we show below, in an irradiated sample 
the Landau quantization leads to violation of the Ein- 
stein relation between the dc conductivity and diffusion 
coefficient. Then a finite photocurrent is driven by a 
built-in electric field even in the sample at a constant 
electrochemical potential. In an open circuit, a photo- 
voltage is produced. Both these photo-galvanic signals 
oscillate around zero as a function of magnetic field as 
observed in the experiment. Another motivation for the 
present study is the physics of ZRS where the uniform 
charge and field distributions become electrically unsta- 
ble, and the knowledge of the transport properties of in- 
homogeneous system is of central importance for determi- 
nation of the configuration and dynamics of the current 
domainsi^i^i^i^ 

The paper is organized as follows. In next section we 
formulate an approach applicable for description of the 
electron kinetics in the presence of non uniformly varying 
potentials. In Sec. |TTT] we discuss the steady state dis- 
tributions and current in the absence of the microwave 
illumination for different experimental setups. In Sec. II VI 
the microwave-induced magnetooscillations in the local 
transport coefficients are calculated. In Sec. |V] we es- 
tablish the relation between the local transport coeffi- 
cients and the photocurrent or photovoltage oscillations 
observed in the experiments. Main findings are summa- 
rized in Sec. IVII 
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II. ELECTRON KINETICS IN COORDINATE 
AND ENERGY SPACE 

We consider a 2DEG in a classically strong magnetic field 
(wcTtr ^ 1, where rtr is the transport scattering time), 
and in high Landau levels (chemical potential /i ^ ujc)- 
Hereafter we put h = 1. Transport of electrons in such 
system is most conveniently formulated in terms of mi- 
gration of the guiding center R(t) of the cyclotron orbits. 
The dissipative component of the dc current is given by 
the rate of changes of R(i) due to collisions with impu- 
rities 



en(5fR)coU : 



(1) 



where n is the 2D electron density, and e = — |e| is 
the electron charge. We describe these collisions using 
a generic disorder model characterized by an arbitrary 
dependence of the elastic scattering rate 



n— — CK 
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(2) 



on the momentum scattering angle ip2 — fi- Every scat- 
tering event is accompanied by the shift of the guiding 
center by AR12 — Rc&z x (ni — n2), where Rc = vpj^c 
is the cyclotron radius, vf the Fermi velocity, and rifc = 
(cos <y9fc, sin (y9fc) the unit vector in the direction of motion. 

Assuming weak one-dimensional spatial variations of 
n — n(x) and of the electrostatic potential </> = 0(a;), we 
express the migration of the guiding center in terms of 
the local distribution function /^j. and the local density 
of states (DOS) v^x = v\e — e^(x)\ where vq = mj^Tx 
is the DOS at i? = 0. In equilibrium, the distribution 
function / depends only on the total energy of electron e 
and is characterized by the position independent electro- 
chemical potential 7y(x) = \x{x) -I- e<l){x) = const(x)i^ By 
contrast, the DOS in high Landau levels is a periodic 
function of the kinetic e nergy s — edjx). Generalizing 
the approach of Refs. [Ulili^liiligllsa to the present 
spatially inhomogeneous case, we obtain 



scattering off disorder are given by 

1 Pv+<p' . 2V -V' 



T.(0l) 
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(6) 
(7) 



The microwave field (screened by the 2D electrons^) is 
taken in the form 



(8) 



where 2^/^e± = ii^y and the complex vector (s-|_, s-_) 
of unit length characterizes the polarization. The dimen- 
sionless power Pg is 

Pg = r - 2Re[8+£*_e'^], (9) 
r = \£+\' + \£^\\ (10) 
£± = s±evpE^uj^^{u} ± ujc)^^ . (11) 

Apart from the modification of the scattering integral, 
the microwave illumination leads to a nonequilibrium en- 
ergy distribution fg^ of electrons, which is controlled by 
inelastic relaxation. The corresponding balance equation 
reads 



fs: 



AT) 
J ex 



(12) 

where Tin ^ T ^/i is the energy relaxation time due to 
electron-electron interaction,^^ Ai is defined in Eq. ([5]), 
x' — X + AX^^,, and/(^) is an equilibrium distribution 
function. As we will see below, Eqs. ^ and describe 
both the displacement and the inelastic contributions to 
photovoltage (or photocurrent) oscillations 



III. DARK STEADY STATE 



X 00 

j:E = 2i^oe J dxi J dx2iWxi^x2 -Wx2~rxi ), (3) 
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Wx,-,X2 = ( Ml\%5{xi -X2 + AX^,^,) 

X {rl%6ie^ ~ 82) + r(ph^)^ '5(^1 - ^2 ± ^)}).(4) 



Here AX^ 



ex • AR12 = i?c (sin — sin (^2) is the x- 



component of the guiding center shift, factor 2 accounts 
for the spin degree of freedom, the angular brackets de- 
note averaging over angles ipi^2 and integrations over £1,2, 
and 



■^1%' = i'exi'e'x'fexil - fe'x'] ■ 



(5) 



The rates of elastic [F^'^'^j] and photon-assisted [F^'^j'^^^] 



A. Infinite 2DEG in a constant electric field 

Equat ions similar to Eqs. ([3]), (|12p were used in Refs. l45l. 
kdkTllislligllsol for analysis of MIRO in a homogeneous 
case of a constant electric field E in an infinite 2DEG. 
In this case, the dark (nonequilibrium due to dc current) 
distribution 



(T) 



exp ■ 



ri{x) 



T 



r]{x) — fj.{x) + e<f){x) 



(13) 
(14) 



is characterized by a coordinate-independent local chem- 
ical potential /i = const (x). The occupation of all states 
having equal kinetic energy e — e(j){x) = e + eEx is the 
same, see Fig. [T^. 
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FIG. 1: Illustration of the steady state distribution f^x 
(greyscale plot) and position of the Landau levels at £jv (x) — 
e<j>{x) + {N + l/2)u>c (dashed lines) in (a) a non-equilibrium 
state with a constant electric field E in an infinite 2DEG and 
of (b) an equilibrium finite 2DEG with a built-in electric field 
near a contact. 



B. Inhomogeneous equilibrium state 

Let us now consider an equilibrium 2DEG in the absence 
of the microwave field but in the presence of a built-in 
static electric field. In the absence of external voltage 
applied to the sample, this static field can be created, for 
instance, by a metallic contact, as in Fig.[T]3. The distri- 
bution function (|13p in this inhomogeneous equilibrium 
case is characterized by a position-independent electro- 
chemical potential 



T] = const(a;). 



(15) 



while both electron concentration and electrical poten- 
tial vary with x. By contrast to the previous case of a 
constant electric field, here all states having equal total 
energy e are equally occupied, /iJ'' = fe'^\ 

In the absence of the microwave field, — 0, the in- 
elastic part oc T'i^^J^ = of the total current ^ is absent 
while the elastic part proportional to Al^/^ — %' — 
vanishes as f^x' = fex, see Eq. ([5]). We arrived at the 
trivial result that the current does not't flow when the 
2DEG is in the equilibrium state characterized by a given 
temperature T and electrochemical potential rj. 



Electrical current and diffusion; Einstein 
relation 



Vanishing of electrical current in an equilibrium state 
of inhomogeneous system can be equivalently formulated 
as the Einstein relation between the linear-response con- 
ductivity and diffusion coefficient. It is instructive to 
derive this relation by considering weak perturbations of 



a spatially homogeneous equilibrium system with a fixed 
concentration of electrons n — const(x) [and, therefore, 
fi = const (x)]. 

According to Eq. ([3]), the electric current induced in 
this system by the infinitesimally small electric field E = 
—V4> reads 



(16) 
(17) 



Here au = e^j/o-Rc/2Ttr is the classical Drude conductiv- 
ity per spin orientation in a strong B (wc^tr 3> 1), the 
transport relaxation time Ttr is expressed in terms of the 
moments r„, Eq. ([2]), as Ttr = {tq^ ~Ti^)~^j and the su- 
perscript "(dark)" refers to the equilibrium state in the 
absence of microwaves. 

Now we put E ~ —V(j) — and calculate the diffusion 
current, i.e. the linear response to a small gradient V?i 
of the concentration 

n{x) = 2uo j deu{e)fex ■ (18) 

The diffusion current 

-eD^'^^'^Vn, (19) 
defines the dark diffusion coefficient 

^(dark) ^ 2l'oX'('*"''^ /X^''"''^ ■ (20) 

Using n given by Eq. (jlSp , we express the dark compress- 
ibility as 

^(dark) ^ Q^^g^^ ^ ^^i^o [ dei>{e)dj^^l (21) 



■(dark) 



■(dark) | 
J '|V0=O 



The quantity pf'^'*'^'^) has the dimensionality of the diffu- 
sion coefficient Z?''^'*'^'^^ and is defined through the relation 

V0=o - -2ej.o2?^''"'''^V/i(x) (22) 



as the current response to the gradient of the chemical 
potential V/i(a;). Calculation using Eq. ^ gives 



p(dark) fdei>^Q^fiT) 
2Ttr J 



(23) 



If we now allow for a generic weak perturbation of a 
homogeneous equilibrium system, the local current takes 
the form 

^■(dark) ^ -2a^'^'''^^\7(l){x) - eD^'^'''^^Vn{x) . (24) 



According to Eqs. (|T6|) - ((23)) . the diffusion coefficient 
£)(dark) conductivity 2a^'^^''^^ are related as 



2^(dark) _ g2^(dark)£)(dark) 



(25) 



The Einstein relation (p5| ensures the absence of the elec- 
tron fiow in the equilibrium state (jl5p with 

Vry = + eV0 = (x''^"''^)"^Vn -I- eVcj) = 0. (26) 
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Related to Eq. (|25|) is a simpler identity 

^(dark) ^ g2^^p(dark) ^ (27) 

which does not involve the compressibility. 

An important consequence of the Einstein relation is 
that the current response of any equilibrium system can 
be represented as 

^.(dark) ^ -^2eiyoV<-'^'"'^Wr]{x) , (28) 

i.e. the current is proportional to the gradient of the 
electrochemical potential independently of what kind of 
perturbation causes the current flow. 

In what follows we assume an experimentally relevant 
range of high temperatures, 2tt'^T/lUc ^ 1, where Shub- 
nikov - de Haas oscillations are thermally suppressed and 
transport properties are independent of the position of 
the chemical potential with respect to Landau levels 
In this limit, Eq. (fTT]) reduces to 

^(dark) ^ ^^^~2^^ ^ 2tT^T/cJ, » 1. (29) 

Here (. . implies energy averaging over the Wc-periodic 
DOS oscillations. In high-T limit, the dark compressibil- 
ity ((2T|) reduces to 21/0, 

x'^"""^ = 2iyo{i>e)e = 2j^o , 2tt^T/lj, > 1. (30) 

A simple linear relation n(x) = 2i'Qfi{x) makes two defi- 
nitions of the diffusion current p9)) and ((22|) identical, 

^(dark) ^ p(dark) ^ i^(^2^^ ^ ^ ^ ^ (3^) 

2Ttr UJc 

In the presence of microwaves, however, the two 
definitions are not fully equivalent in view of the 
microwave-induced magnetooscillations in the compress- 
ibility (MICO),''^ see discussion in Sec. IIVCI and in 
Sec. ESI 

IV. LOCAL CONDUCTIVITY AND DIFFUSION 
COEFFICIENT IN ILLUMINATED 2DEG 



A. Nonequilibrium current flow 

We now turn to evaluation of the transport properties in 
the presence of microwave radiation. The key observation 
is that in the nonequilibrium steady state the Einstein 
relation (P7)l between the dc conductivity a and diffusion 
coefficient T) does not hold anymore, a ^ e^v^D. In other 
words, the current cannot be represented in the form of 
Eq. (pS)) with some modified transport coefficient and 
electrochemical potential. According to our calculation 
based on Eqs. ^ and p2p . the nonequilibrium dc current 

= -2cr^V0 - 2evQVSI^. (32) 



necessarily contains an extra "anomalous term" — 2cr.,,V(/> 
violating the Einstein law. In these terms, the total 
conductivity cr, which defines the dc current j = 2a E 
in a homogeneous system (as in the case of MIRO, see 
Sec, nil A| . is given by 

(7 = cr,, 4- e^i'ol?, (33) 

while the diffusion coefficient D entering the current 
jIb^o = —eDVn at i? = 0, is expressed through the 
nonequilibrium compressibility Xj~ 

D = 2yaX'^V, (34) 

similar to Eq. ((20|) . In next two subsections we calculate 
the anomalous conductivity ct^ and a photoinduced part 
of the diffusion coefficient to the minimal order E"^ . 

B. Anomalous component of conductivity 

In this subsection we calculate the anomalous component 
of the conductivity. For that purpose we put Vry = 
and use Eqs. ([3]) and (fT2|) with the position-independent 
dark distribution, 

fP = [e^^-")/^ + r? = const(x). (35) 

Similar to the dark case, Sec lIIIB] the microwave correc- 
tion ^ to the elastic scattering rate gives no contribution 
to the current ([3|) due to cancellation M%,% — M^, ^, = 0. 
Therefore, the current cx E"^ can be either due to (i) the 
microwave-assisted scattering off disorder (represented 
by the second term oc V^'^^^ in Eq. ^ with unperturbed 
f — fe , displacement mechanism) or due to (ii) the po- 
sition dependence of the microwave-induced nonequilib- 
rium distribution ()12p [modifying the clastic term in the 
total current ^ , inelastic mechanism] . Correspondingly, 
the anomalous conductivity is a sum of the displace- 
ment and inelastic contributions 

a, = -f ) (36) 

Displacement contribution (t^'^'*'\ Using the position- 
independent distribution ([35]), the second term of Eq. ([3]) 
can be represented as 

4-) = 2e.oE / deifP - /iE)(e(A)rSJ? 

X 

X J dx' v[e - e(/)(x')]j>[£ ± w - e(/i(x' + A)])^^/(37) 

where the Heaviside function 6(A) imposes the condition 
A = l^X^^p- > 0, see Eq. The electric field enters 
this expression only through the position dependence of 
the local DOS i>i;x — v[e — e(j){x)]. In the absence of the 
local electric field E = —V4>{x), two terms of Eq. ([ST)) 
corresponding to e ± w exactly cancel each other. The 
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terms linear in E produce the displacement contribution 
to the current j'^^^'^ = a^^'^^^^ E under condition (fT5|) : 



Microwave-induced oscillations of the diffusion 
coefficient 



fif = aD-^(P - Re[£+£*])7^l(a;), (38) 



Here we used I-k^T juj^ 1, and 



r-i=3ro-i-4rri+r2-i 



(39) 



(40) 



in terms of Eq. Function 7?,i(tj) oscillating with 

the ratio iojujc is specified in Sec. IIVDI for two limits 
of strongly overlapping and well separated Landau lev- 
els, together with similar oscillatory functions entering 
Eqs. dm) and ([iTjl . 



Inelastic contribution cr^'"'' . The inelastic contribution 
is due to microwave-induced changes Sf in the distribu- 
tion function. To the leading order E'^ and in the limit 
V0(a;) -> 0, Eqs. ^ and ^ give 

Sfe. = ^(/f ) - fgl)i>[e ±u;- e^x)]. (41) 



In contrast to the spatially independent /i^"*, the 
microwave-induced part Sf of the electronic distribution 
oscillates in coordinate space at a fixed total energy e due 
to spatial oscillations of DOS i>ex in Landau levels tilted 
by the electric field. As a result, the elastic contribution 
to the current ^ does not vanish, M%,% — M% %, ^ 0. 
Substitution of Eq. (|4T|) for f^^ in the elastic term of 
Eq. m with rj;l : 
tion to the current. 



T^l^' produces inelastic contribu- 



,-(in) _ 



^2euo I de j dx' ^'-'^^'^^^ 

J \ J Tip — ipf 



(42) 



where A = AX^p^p/ > as above in Eq. ([57]) . Assuming 
I'K^Tjuic 3> 1 and keeping the linear term oc V</)(x), we 
obtain j,'/"^ = a'r\°'^ E with 



4rtr 



(43) 
(44) 



It is worth mentioning that both a,',*^'^^ and cr,^/"^ orig- 
inate from the spatial dependence of the DOS which re- 
quires both the Landau quantization and the presence 
of electric field. In the absence of Landau quantization, 
i> = 1, functions TZi{ll!) and 7^2(ijj) entering Eqs. ([39| and 
vanish (see Sec. lIVDl '). Therefore, within our model 
the Einstein relation of the dc conductivity and diffusion 
coefficient is restored in the classical^ limit uJcTq 0: 



and 



Jin) 



0, see Eqs. 



Now we assume Wcj) = and Vfi — Vrj ^ 0, and calculate 
the microwave-induced correction Dph = D — Z?''^'*'''^^ to 
the diffusion coefficient, see Eqs. (p2)) and ((34| . In con- 
trast to the previous subsection, now the DOS is po- 
sition independent, while the dark distribution varies 

(T) 

in space. Using the linear approximation fl^\^sx ~ 

feV —SxV^defex^ in Eq. ([3]), we obtain jj; — — 2et'oPV/i 
with 

j de[-d,fiV]{Q{AX^^,)AXl^, 

X [ulTf;i,+i.A^,^^ + ^e+u)T^^$])^^'. (45) 

Performing the angular and thermal averaging for 
2tt'^T/uJc > 1, we get, similar to Eqs. 

2?ph = ^{V- M£+£'L])ni{uj), (46) 



(47) 



7^3(w) = {vl - t>eZ>e+^)e , 



where Dph = I?-r>(dark) ^nd p(dark) jg gjygj^ by Eq. jST]). 

In Eqs. ((45|) and (j46ll . the microwave- induced changes 
in the distribution function (|4ip were not taken into ac- 
count. The reason is that in the limit 27r^r/wc ^ 1 
the corresponding contribution to the diffusion coefficient 
is exponentially suppressed. The inelastic contribution 
obtained from Eqs. (gl]) and ^ using V0 = 
and V/i ^ 0, reads 



V 



(in) 
ph 



devlv,±^d,\fiV-lilix]- (48) 



In the limit ^tP'T juic ^ 1, this expression vanishes simi- 
lar to the Shubnikov-de Haas oscillations. Therefore, the 
microwave-induced oscillations in the distribution func- 
tion (|41]) produce the contribution (|43|) to the anomalous 
conductivity ct^ only, while the displacement mechanism 
provide similar oscillations both in ci,;, Eq. (j38p . and in 
P, Eq. (gSl). 

While the dark quantities p(daik) ^^^^ j:,(d^rk) ^^^^ 
identical at 2n'^T/uJc ^ 1, see Eq. ([51]) . in the pres- 
ence of microwaves they are not equivalent in view 
of the microwave-induced compressibility oscillations 
(MIC0).'^2 MICO do not enter the quantity V since it 
is defined through V/i, but modify the diffusion coeffi- 
cient D defined through Vrt = x^/^i which gives 



D 



ph 



D-D 



(dark) 



ph 



(49) 



However, since we are interested in linear-in-i?^ correc- 
tions to and in the present work, we 
can approximate the compressibility by its dark value 
^(dark) ^ 2jyo (thus neglecting MICO that lead to terms 
oc E* in -Dph). Moreover, even at high orders in E^, the 
compressibility can be approximated as x = dn/dfj, = 



6 



2fo assuming spatial variations of n{x) are smooth on 
a scale of the inelastic length. At shorter length scales, 
MICO can be strong (of order i/q). This situation arises, 
in particular, in the regime of ZRS,^^ see Sec. IV Al 

So far we considered the two cases Vrj — and V(j) — 
which give, correspondingly, the anomalous conductivity 

and the photoinduced correction to the diffusion coef- 
ficient D = D^'^'"'^^ + Dph. The sum cTph = a^ + e'^voVpi,, 
given by Eqs. (gS]), and reproduces 

the result o'^^'^^i"^'''''^^i'^^'^° obtained earlier for the homo- 
geneous case of the MIROj^^'^'^'^i^ which corresponds to 
V/i = and to the constant — V0 = —Vrj/e = E, see 
also Sec. IVC] 



D. Form of the oscillations for overlapping and 
separated Landau levels 

The form and the phase of the magnetooscillations in the 
anomalous conductivity (t,, and in the diffusion coefficient 
D, Eqs. ([39]), dB]), and ([it]), as well as the quantum cor- 
rection to the dark conductivity, Eqs. are expressed 
through the certain energy averages 7?.„(cj) over the pe- 
riod ojr of the DOS oscillations, 



^1 II"g"(^")^2(|^^«|), 



/.~.2/~. 



TZi{u!) = ujd^ 

7^.2 (w) = ajdu^iv^ive+u, + v, 
7^3(w) = {vl - 



))e 



(50) 



Here we specify these functions in two limits of strongly 
overlapping and of well-separated Landau levels (LLs) 
within the self-consistent Born approximation (SCBA). 
At high LLs, sp 3> ti^,Wcj disorder can be treated within 
the SCBA provided the disorder correlation length sat- 
isfies d <^ Ib and d <C vpTq, where Ib — (c/ei?)^/^ is 
the magnetic length and Tq the quantum relaxation time 
[rq = To in terms of the moments r„, Eq. , 45,46,47,74 

In moderate magnetic field, 6 = exp(— 7r/a;crq) <C 1, 
LLs strongly overlap and the DOS is only weakly mod- 
ulated by the magnetic field, i>{e) = 1 — 2(5cos(27re/a;c)- 
In this limit 



7^o = 1 + 2S^ 



7^l = -2S' 



, 2ttuj . 27rti; 



■ sm ■ 



UJc 

, 2i:lj , 2t:uj 



■ sm ■ 



Us = 4(5^ sin^ 



UJc UJc 

nuj 



(51) 



In the limit of separated LLs, uJcTq ^ 1, the DOS is a 
sequence of semicircles of width 2r = 2(2u;c/7'"''q)^^^ 
UJc, i.e., D{e) — TqRe-^/F^ — (&)2, where 6e is the de- 
tuning from the center of the nearest LL. In this limit, 
calculation yields 



AuJUJr X — > 

7^2 - -7^o^;5-2]sgn(^!„)$2(|^^„ 

n 

7e,, = 7^n i-^Hidaj) 



(53) 
(54) 

(55) 



7^0 = WUc/iTT^T, 



(52) 



The parameterless functions of ri„ = (uj — nuJc)/T are 
nonzero at 0< |r2„| <2, where they are expressed as 

Hiix) = i2 + x)[{A + x^)E{X) -AxK{X)]/8, (56) 
H2(x) = 3a; [ (2 + x)E{X) - 4K{X)]/8 , (57) 
Att^2{x) = 3a;arccos(a; - 1) - x{l + x)^ x{2 - x) . (58) 

Here X = {2 — xYl (2-1- a;)^ and the functions E and K are 
the complete elliptic integrals of the first and second kind, 
respectively. Graphical representation of the functions 
(IMll-dMl) can be found in Ref. [tBI. 

In the crossover magnetic field, ujcT^ ^ 1, functions 
([50)1 obtained using analytical expressions for the DOS 
become very cumbersome. In this crossover region, the 
form of the oscillations can be obtained using numerical 
solution of the SCBA equations.^" In particular, such nu- 
merical solution is used in the calculation illustrated in 
Fig. ^Bp below. 



V. PHOTOCURRENT AND PHOTOVOLTAGE 
OSCILLATIONS 

In this section, we use the obtained local transport co- 
efficients for calculation of the electrical current in dif- 
ferent experimental situations. The effects related to 
microwave-induced modifications of the spatial distribu- 
tion of carriers and fields are considered in Sec. IV Al 
The current-voltage characteristics (CVC) of an infinite 
2DEG stripe between two metallic contacts are obtained 
in Sec. IVll Using these CVC, in Sec. |VS] we calcu- 
late the photocurrent and photovoltage and compare our 
findings with the experiment of Ref. [gJ. In Sec. IV Dl 
nonlinear effects in the photovoltage (with respect to the 
microwave power) are discussed, which were observed 
experimentally^ and are also well reproduced by the the- 
ory. 



Photo-induced changes in the field and charge 
distribution 



In the presence of nonuniform carrier and field distri- 
butions, as in Fig. [T|d, the local transport coefficients 
an and V entering the local current ([5^ do not com- 
pletely determine the transport. The full theory should 
include a self-consistent solution of the Poisson and con- 
tinuity equations for a given experimental setup. Indeed, 
the photoinduced current density, given by Eq. (j32p with 
the dark profile of the electrostatic potential (/)(^^'''^) and 
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with Vrj = 0, is jx = — 2(t^V(/)'-'^^'''^-'. In general, such 
jx does not satisfy the continuity equation Vj^ = 
in view of a nonlinear spatial variation of 0^'^"''^ as, 
for instance, in Fig. [T|d. Therefore, the photoinduced 
variation of the electron density 5n{x) = n — ridark (we 
assume |(5n(x)| ^ n) and of the electrostatic potential 
6(f>{x) = 4> — (/)('^^'''^) should be taken into account. The 
latter are related to each other by the inverse capacitance 
matrix W as 

5<Px^ejd.Wx.8nx. (59) 
Using the relation 

Srj ~ e5(j) + x^^Sn (60) 

valid at 2tt^T /lUc ^ 1, we represent the Poisson equation 
in the form 

5ri=[l + {e^xW)-^]e5(j). (61) 

Using Eqs. (l32|) and (f6T|) for a fixed current density jx — 
j = const(a;), one arrives to a formal solution for the local 
variation dri{x) of the electrochemical potential, 

St] = -e|2CT,,V[l + (e\M?)"i]"i+2e2i/oPv| ^ 
X (j + 2(7^V0('^^^'^)). (62) 

Solution of the above non-local equation is required 
if the amplitude of oscillations in cr^ becomes of or- 
der (T^'*'*'''^) = e^voV'-'^'^''^^ [otherwise one can neglect the 
photoinduced changes Sn{x) in view of the smallness of 
oc E^]. In conventional magnetoresistivity experi- 
ments this corresponds to the regime where the zero- 
resistance states are formed^ '^^i^^i^^i^^i^^ [both theory 
and experiments show that the ZRS appear still in the 
linear regime in the microwave power where Eqs. p6|) . 
([38| , ([43^ and ([46]) still apply] . According to the theory 
of Ref. |43, the ZRS is a manifestation of a spontaneous 
symmetry breaking of a homogeneous state with nega- 
tive resistivity leading to the formation of the current 
domains. In this picture, the residual resistivity in the 
ZRS, which is observed in part of experiments, is due 
to the electron transport across the domain walls and 
near the boundary of the 2DEG. Inside the domains, the 
transport is dissipationless. 

The boundaries of the domains are characterized 
by strongly nonuniform carrier and field distributions. 
Therefore, the results of the present work, in particular, 
the violation of the Einstein relation and the appearance 
of the anomalous component of conductivity ct,,, should 
play an important role for development of microscopic 
theory of transport in the ZRS regime. 

B. Boundary conditions and current-voltage 
characteristics 

We now consider the photocurrent and photovoltage os- 
cillations in a 2DEG with metallic contacts. As we show 



below, specific boundary conditions ([64)) at the interface 
with metallic contacts make the details of the poten- 
tial and carrier distributions inside the sample irrelevant 
[thus, one need not solve a complicated electrostatical 
problem ([62]) ]. More precisely, as long as simple ID or 
Corbino geometry is considered (see Fig. [3]) and the built- 
in electric field is not too strong , the current and voltage 
between the contacts are fully determined by the differ- 
ence of the work functions of the contacts and by the 
local transport coefficients cr^ and V, see Eq. (|65p below. 

Indeed, using the fact that in the linear approximation 
with respect to the dc field not only jx but also and 
V are position independent, and integrating both parts 
of Eq. ((32)1 along a contour connecting two contacts at 
X — and x — we obtain the relation 

jL - 2a^[0(O) - 0(L)] + 2eiyoV[ri{0) - ry(L)]. (63) 

It is natural to assume that microwave radiation does not 
change the electron concentration on the metallic side of 
the interfaces due to a huge density of states there. Since 
both the electrochemical and electrostatic potentials are 
continuous at the interface, this fixes the chemical po- 
tential in the 2DEG near the interfaces. Introducing the 
voltage V = [ri{0) — ri{L)]/e and the difference of the work 
functions of the two contacts eUc^ we write the boundary 
condition in the form 

0(0) - 0(L) - y = [^l{L) - ^l{0)]/e = U,. (64) 

Equations ([55]) and (|64p yield the desired current-voltage 
characteristics (CVC), 

jL^2(j^Uc + 2aV, (65) 

where 2a — 2tT^ -I- 2e^VQ'D is the total conductivity. We 
emphasize that the CVC retains the form (p5)) for ar- 
bitrary microwave power (provided the transport coeffi- 
cients (T,, and V are calculated to all orders in V). Also, 
Eq. (|65p is applicable in the case when the microwave- 
induced redistribution of carriers is significant, an ~ 
^(dark)^ sec Eq. ([5^ [provided the relative change of the 
electron density across the sample remains small]. The 
CVC (|65p is modified only when the linear approxima- 
tion with respect to the dc field breaks down. For such 
strong dc fields, the transport coefficients a^ and D be- 
come field- and coordinate-dependent (and, therefore, are 
no longer uniquely defined). Only such strong dc field 
makes important the details of the electrochemical and 
electrostatic potential distribution in the interior of the 
sample, which necessitates the full solution of the Poisson 
and continuity equations with the boundary conditions 

(El. 

C. Photocurrent and photovoltage 

If the geometry of two contacts is identical and the dif- 
ference of the contact potentials is zero, Uc = 0, the CVC 
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(|65)) reproduces the Ohm law in the bulk, j = 2(tV/L. 
Here a contains the displacement and inelastic contribu- 
tions to the MIRO, Eqs. jHS]), dSH]), gSl), and (gS]), repro- 
ducing the results of previous calculations."*^'— lili^i^i^ 
An asymmetric contact configuration results in a non- 
zero average electric field E ~ Uc/L inside the sample 
in the absence of the bias voltage, V = Q. In the pres- 
ence of the microwave induced anomalous conductivity, 
(Jn 7^ 0, the built-in electric field E = Uc/L leads to the 
photocurrent at zero bias voltage, 



or, in the open circuit, to the photovoltage 



(66) 



(67) 



as observed in the experiment.^'* Two experimental traces 
of the photocurrent for different temperatures are shown 
in Fig. [2^. FigureOj illustrates the inelastic contribution 
(|43| to the anomalous conductivity (|36p. which demon- 
strates an excellent agreement between the theory and 
experiment. A typical small shift of the zeros of the 
photocurrent from the integer and half-integer values of 
uj/ojc in experimental traces is similar to observationsfi>S 
for MIRO and can be attributed^ to a slight deviation 
of the electron effective mass from the standard value 
TO = 0, 067mo used in Fig. [2^. 

In the case of overlapping LLs, the phase and the 
form of the photocurrent oscillations is identical for 
the displacement and inelastic contributions to Gr^, see 
Eqs. dini),®, (HSl), and Therefore, one can distin- 
guish between them only owing to a strong temperature 
dependence of the inelastic scattering rate. The temper- 
ature dependence in Fig. [2^ shows that the inelastic con- 
tribution to the anomalous conductivity cr^ is substan- 
tial. At the same time, this dependence is weaker than 
cj" oc Tin oc predicted by the theory^ at the leading 
order in both the dc and microwave fields, see Eq. 
The weaker T-dependence can be attributed either to a 
strong admixture of the displacement contributio n^^'^^i^° 
at r = l.'bK or to nonlinear effect^ i'^^i^^i'*^ (in the mi- 
crowave power or in the dc field) . Alternatively, it can be 
the manifestation of a noticeable heating of the electron 
gas^^'^'^ at T = O.hK , since the inelastic scattering time 
is a function of the electron temperature rather than the 
bath (phonon) temperature. 

Let us emphasize that the analysis of the ID geome- 
try considered above (Fig. [3^) is directly applicable to 
the case of the Corbino geometry (Fig. ^p) for which 
the Hall conductivity does not enter the local relation 
([5^ . In the case of Corbino geometry, the current den- 
sity j{r) cx r~* is inversely proportional to the distance 
from the center, so that the total current J = 2TTrj{r) is 
conserved. Therefore, integrating both parts of Eq. ([32]) 
along a contour connecting two contacts at r = a and 
r = b one obtains a modified CVC in the form 




0.50 



T = 0.5K 



03 









(b) 











0.0 



0.5 



1.0 



CO / CO 
c 



FIG. 2: (a) Photocurrent between the Corbino-like internal 
and strip-like external contact (see Fig. [S}:) vs. magnetic field 
measured in experimental setup of Ref. [6J for T = 1.5 K 
(solid line) and T — 0.5 K (dashed line); (b) Photocurrent 
[calculated using numerical solution of the SCBA equations 
according to Eqs. p2p and H42|) . see Ref. for details] for 
ujTq = 10 and for a linear polarization of the microwaves. Here 
we took into account the screening of the incoming radiation 
by the 2DEG which results in a strong S-dependence of the 
internal microwave field E^i [entering Eq. (|8])] in the vicinity 
of the cyclotron resonance.— 



J-!- In - = 2(T„ Uc 
2n a 



2aV. 



(68) 



Comparison of Eqs. ([55]) and (|65p shows that results for 
ID geometry (Fig. transform into results for Corbino 
geometry (Fig. [21d) after replacement jL J^^^a- 
In both cases of ID and Corbino geometry, the photo- 
galvanic signals result from non-zero average built-in 
electric field which requires either different work func- 
tions of metallic contacts or different electron densities 
under capacitively coupled gated probes.~ A further pos- 
sible source of the asymmetry is different geometry of the 
contacts, e.g. the Corbino-like geometry of the internal 
contact and the strip-like geometry of the external con- 
tact located on the perimeter of the sample, see Fig. [2t 
(as in the part of the experiment of Ref. |6J that utilized 
heavily doped ohmic contacts). For such geometry, the 
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FIG. 3: Illustration of diflFerent contact geometries, (a) One- 
dimensional (ID) geometry: 2DEG between two long strip- 
like contacts, (b) Corbino geometry: as long as the axial 
symmetry is preserved, can be reduced to the ID case, (c) 
Combined geometry used in the experiment of Ref. [6^ . The 
experiment shows that the photogalvanic signal forms near 
the internal contact having the Corbino geometry while the 
external strip-like contact plays no role. In such conditions, 
the theory for geometry (b) is also applicable to the case (c) . 



photo-galvanic signals were shown^ to be formed in the 
vicinity the internal Corbino-like contact and the above 
consideration should be valid if one puts the difference of 
work functions of the doped internal contact and 2DEG 
instead ofUc in Eqs. §5^, and (p7|) . 



D. Nonlinear effects in the photovoltage; 
Photoresistance 



The magnetooscillations in the photocurrent ([M)) are 
fully determined by the anomalous conductivity cr,, and, 
therefore, the oscillations are symmetric with respect to 
the average value j = 0, see Fig. [21 By contrast, the ex- 
perimental traces of the photovoltage oscillations show^'' 
a strong asymmetry with respect to V^h = value. This 
asymmetry is due to additional microwave-induced oscil- 
lations in the denominator a — o-t^^''^' + + e^foPph 
of Eq. dST]). From previous studies of the MIRO'^^-''* it 
is known that contributions to ct^ of second order in the 
microwave power are still small when the magnitude of 
the first-order terms approaches the dark conductivity 
^(dark)^ This legitimates the use of Eq. ([S7|) in the nonlin- 
ear regime. Neglecting inessential correction 2?ph [which 
is a factor t:lu/uJc smaller than the displacement con- 
to (T„ even in the absence of the inelastic 



tribution 
contribution 



], one can rewrite Eq. (p7)) as 



Vr 



ph 



1 + cr(dark) /^r^ 



(69) 



Equation (|69|) explains a strong asymmetry of the pho- 
tovoltage oscillations observed in the experiment.®^ Fur- 
ther, the nonlinearity of Eq. makes possible the ex- 



perimental determination of the value of the contact po- 
tential difference Uc, since in the formal limit |cr^| 3> 

fj(dark) j^g^g simply Vph = —Uc- 

Apart from the photocurrent and photovoltage, one 
can measure the two-point differential photoresistance 
Rph — dV/dl by driving a small current through the 
sample. Such measurements were also done in Ref. Is^ . 
and the results were compared to the ratio Vph//ph taken 
from two independent measurements of Vp^ and /p^. The 
comparison demonstrated very good agreement, as can 
be expected from Eq. ((65)) giving 



dV 
'dj 



L 
2^ 



and Eqs. and ([S7)) yielding 

jph 2cr 



(70) 



(71) 



The measured photoresistance showed clear magne- 
tooscillations with the phase opposite to the MIRO. The 
phase shift of oscillations by tt is in agreement with 
Eqs. ((70| . (|7T|) predicting i?ph oc cr~^, which should be 
compared with p^x oc in conventional magnetoresistiv- 
ity measurements of the MIRO J^iS^^^^^ 



VI. CONCLUSION 

Summarizing, we have presented a quantum transport 
theory for a 2DEG in high Landau levels illuminated by 
the microwave radiation in the presence of a spatially in- 
homogeneous dc electric field. The theory explains the 
microwave-induced photocurrent and photovoltage oscil- 
lations observed in the recent experiment.— 

We have shown that in an irradiated sample the Lan- 
dau quantization leads to violation of the Einstein rela- 
tion between the dc conductivity and the diffusion co- 
efficient. As a result, a non-zero average electric field 
leads to the electric current which is not compensated by 
the diffusion flow even for the electrochemical potential 
remaining constant in space. The experimental observa- 
tion of the effect requires an asyinmetry (for instance, 
in contact geometry, as in Ref. |64| . or in material com- 
position of two contacts) which determines the direction 
of the current. At the same time, the obtained current- 
voltage characteristics are shown to be independent of 
detailed potential profile in the sample provided the rel- 
ative change in the electron density across the sample 
remains small. 

The effects discussed in this work should also play 
an essential role for the transport in the zero resistance 



states 



3.4.5.6 



In this regime, the uniform charge and field 
distributions become electrically unstable4^ The sys- 
tem breaks into current domains and peculiarities of the 
transport properties of the inhomogeneous system be- 
come of central importance. 
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